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Abstract. We study an interacting electronic system exhibiting a spin nematic 
instability. Using a phenomenological form for the spin fluctuation spectrum near 
the spin-density-wave (SDW) phase, we compute the spin nematic susceptibility in 
energy and momentum space as a function of temperature and the magnetic correlation 
length £. The spin nematic instability occurs when £ reaches a critical value £ cr , i.e., its 
transition temperature Isn is always higher than the SDW critical temperature Tsdw- 
In particular, £ cr decreases monotonically with increasing Xsn- Concomitantly, low- 
energy nematic fluctuations are present in a wider temperature region as Tsn becomes 
higher. Approaching the spin nematic instability, the nematic spectral function at zero 
momentum exhibits a central peak as a function of energy for a finite temperature and 
a soft mode at zero temperature. These properties originate from the general feature 
that the imaginary part of the spin-fluctuation bubble has a term linear in energy and 
its coefficient is proportional to the square of temperature. Furthermore we find that 
the nematic spectral function exhibits a diffusive peak around zero momentum and 
zero energy without clear dispersive features. A possible phase diagram for the spin 
nematic and SDW transitions is also discussed. 
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1. Introduction 

The mechanism of high-temperature superconductivity is one of the major issues in 
condensed matter physics. By carrier doping into antiferromagnetic Mott insulators, 
cuprate superconductors attain high critical temperatures T c , typically of more than 77 
K, the boiling temperature of liquid nitrogen at ambient pressure. Since it is widely 
accepted that doped cuprates can be described within a one-band model it is natural 
to assume that the important fluctuations are spin fluctuations in these systems. In 
2008 another family of high-T c superconductors was discovered consisting of iron-based 
pnictides [T] . Similar to the cuprate case, superconductivity occurs close to a magnetic 
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phase. However, in contrast to the cuprates, the magnetic phase is metallic and several 
bands cross the Fermi energy, suggesting the importance of orbital degrees of freedom. 
Thus there are two different kinds of fluctuations in the pnictides which may be relevant 
for their physical properties, namely, spin and orbital fluctuations. 

In general it is not easy to assess the relative importance of the two kinds of 
fluctuations in observable quantities. It is known that superconducting order parameters 
due to spin fluctuations tend to have s± symmetry whereas those due to orbital 
fluctuations have s ++ symmetry mn- While this difference in the gap symmetries 
can contribute to identifying the underlying mechanism of superconductivity, a recent 
theoretical study [ 6 ] shows that s± pairing gap can be stabilized even for orbital 
fluctuations when a partial contribution from spin fluctuations is taken into account. 
Experimentally the superconducting phase occurs often closer to the nematic than to 
the magnetic phase mmm- From this one may conclude that nematic fluctuations 
play also an important role for superconductivity and more generally for the physics of 
pnictides. It is thus important to study the properties of nematic fluctuations and their 
origin in detail. 

The nematic instability is accompanied by a structural phase transition from a 
tetragonal to an orthorhombic phase. This structural phase transition is believed to 
be driven by the coupling to electronic degrees of freedom | 8 j, because the observed 
anisotropy of resistivity [ 10 ], the optical conductivity mi and the splitting of electronic 
bands [12] is much larger than what one would expect from the lattice anisotropy. Due 
to these effects the orthorhombic phase in the pnictides, which breaks the C 4 symmetry 
of the normal state, is referred to as an electronic nematic phase. Depending on the 
electronic degrees of freedom responsible for the nematicity, three kinds of nematicity 
can be distinguished: charge da na nans], orbital PH HU, and spin nematicity [19] . 
For iron-based pnictides, the latter two possibilities seem to be relevant. 

Orbital nematicity in pnictides is associated with a spontaneous difference in the 
occupation of d yz and d zx orbitals [201 ED EEJ- Below the transition temperature the 
orientational symmetry is broken, whereas all the other symmetries remain unbroken. 
Orbital fluctuations may lead to a superconducting state driven mainly by fluctuations 
around zero momentum [4]. In this scenario, nesting of the Fermi surfaces is not crucial 
for pairing. Orbital fluctuations around a finite momentum PI E] may, similarly 
to spin fluctuations m si, also lead to superconductivity but in this case nesting 
properties of the Fermi surface are important. Considering general orbital fluctuations, 
retardation effects and quasi-particle weights within Eliashberg theory it was shown 
[24] that orbital nematic fluctuations produce strong coupling superconductivity. The 
obtained transition temperatures and the fact that superconductivity also occurs inside 
the nematic phase are compatible with the experimental observations. 

Spin nematicity in iron-based pnictides is associated with a spontaneous anisotropy 
of the spin fluctuation spectrum in momentum space between ( 7 r, 0) and (0, 7 r) [25j [26] . 
Spin nematicity owes its existence to spin fluctuations and is caused by the interaction 
between them as reviewed in [27] and [28]. It has the general property that it can occur 
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above the spin-density-wave (SDW) phase |27j . which is in nice agreement with the 
experimental findings. 

In general it is difficult to distinguish between spin and orbital nematicity because 
there is a coupling between the corresponding order parameters. If the orbital (spin) 
nematic instability occurs, the spin (orbital) nematic order is induced via a linear 
coupling between them [28] . Hence both orbital and spin nematic orders occur always 
at the same time. Nevertheless, one may argue that the original instability occurs only 
in the orbital or the spin section. As a result one mechanism will dominate and it makes 
sense to study the two kinds of mechanisms leading to the nematic state independently. 

Using linear response theory we will investigate in this paper properties near the 
spin nematic instability. The results are rather independent of microscopic details 
because the functional form for the spin fluctuation propagator becomes generic close to 
the SDW phase. The characteristic features of the low-energy spin nematic fluctuations 
are determined by the anomalous asymptotic behavior of the spin fluctuation bubble 
at low frequencies. Our results for the dynamic structure factor are predictions which 
can be checked by inelastic light scattering. We also present a schematic phase diagram 
containing spin nematic and magnetic transition temperatures. 


2. Model and formalism 


The spin nematicity in iron-based superconductors can be formulated in terms of an 
action as in [29]. Here we wish to formulate it in terms of a conventional operator 
formalism and to clarify the diagrammatic structure of spin nematic physics. 

Iron-based superconductors are often described by a five-band Hubbard model [3|. 
To describe the spin nematic interaction in a microscopic model, we focus on the effective 
interaction of the total spin operator. Then our microscopic model consists of electrons 
hopping between the sites of a square lattice and interacting via their total spins. The 
corresponding Hamiltonian becomes 

H =Y1 e ka C k,a,a C k,a,a + ^ J ( C l) S (<l) ' S (-<l)> W 

k,a:,cr q 

with the total spin operator 


S(q) 


_ t 

9 / , C k,a,a T w lC k+q,a,a'- 


k,a,cr ,cr' 


( 2 ) 


c kaa creates an electron in the band a with spin direction a, and Ck,a, CT annihilates this 
electron, r is the vector of the three Pauli matrices, ek a the one-particle energies of the 
band a, J(q) the strength of the spin-spin interaction in momentum space, and N is 
the total number of the lattice sites. 

Besides the total spin operator we will consider the spin nematic operator 

$ ( q ) = 7fiv^ 7k s ( k + q / 2 )' S(_k + ( i/ 2 )- 


(3) 
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7k is equal to one if k is close to (7r, 0 ) or (—77 0) and equal to minus one if k is close to 
( 0 ,7T) or ( 0 , —7r). 7k has d- wave symmetry and the sum over k in Eq. ([ 3 ]) runs over the 
entire Brillouin zone (BZ). 

The dynamic spin susceptibility is defined by 

i r°° 

x(q ,u}) = j^j ^([SKq^),S z (-q,0)])e^, (4) 

where the index l denotes a Cartesian component, [•, •] the commutator, and 00 includes 
tacitly a small imaginary part 77; we consider a state with spin rotational symmetry and 
thus y does not depend on l. I11 a similar way we define a spin nematic susceptibility 
by 


j r°° 

Xsx(q ,u) = —J df([$(q,f),$(-q,0)])e^. ( 5 ) 

Going over to the Matsubara representation, the usual diagrammatic perturbation 
expansion for electronic systems holds for the above two correlation functions. 
Concerning y(q, u) we are interested in the following in its low-frequency and long- 
wavelength behavior. Here long-wavelength means momenta near (77 0 ) or ( 0 ,7r) or 
the equivalent points in the BZ. Thus we do not try to calculate y by a perturbation 
expansion but expand the inverse susceptibility in powers of to and q — Q q where Q q 
denotes (7r, 0 ) for q close to (7r, 0 ) and ( 0 , 7 r) for q close to ( 0 , 7 r). The result is 

_c/7_ 

r + (q-Q q ) 2 -707/7' 

r is equal to l /£ 2 where £ is the correlation length. 7 is a damping constant and c 
determines the spectral weight. The above parametrization of y is a general form of 
spin fluctuations near the SDW phase and in fact describes rather well the measured 
imaginary part of the spin susceptibility in BaFei.85Coo.i5As2 by choosing appropriate 
parameters p 30 j. 

The lowest order diagrams for the evaluation of ysN are shown in figure [Q(a). The 
wavy line stands for the spin propagator —y, i.e., for a two-particle propagator. The 
first diagram represents y^j/ for the non-interacting case, whereas the second term comes 
from the spin nematic interaction. The effective interaction (solid square) represents 
the four-spin vertex as shown in figures Q](b) and (c). The open square denotes a ring 
diagram where the sides of the square represent electronic Green’s functions and the 
corner points of the square the interaction term in Eq. ([T]). Taking the low-frequency, 
long-wavelength limit in external spin legs the four-spin vertex collapses to one constant 
g. The effective interaction term assumes then the form 

= i)4>(-q)- (?) 

q 

The resulting diagrams for ygN are shown in Fig. 1 . Altogether one obtains 



xsN(q,w) 


xsN(q^) 

1 -f/Xs°N(q> w ) 


( 8 ) 
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Figure 1. Graphical representation of the spin nematic susceptibility \sn- The wavy 
line indicates spin fluctuations and the vertex with a circle in (a) represents a form 
factor 7 k associated with the spin nematic order parameter [see mi The solid square 
is an effective four-point vertex of spin fluctuations defined in (b). The open square 
in (b) and (c) is a bare four-point vertex of spin fluctuations due to the spin nematic 
interaction [see ©]• 


The analytic expression for Xsx is 

Xs°N(q, iu n ) = — ^ x(k + q/2, iu n + ii/ m )x(- k + q/ 2 , ~^m), (9) 

k,m 

where c o n and u m are bosonic Matsubara frequencies, the factor 3 comes from spin 
rotational symmetry and T is the temperature. In the long-wavelength and low- 
frequency limit follows from the above equation for small r, 

Xsn(°. °) = jE W k ’ iVm ^ 2 [(x(k °)] 2 ~ 7' ( 10 ) 

k.m k 


Assuming g to be positive [attractive interaction in (J7|)], (ITU]) implies that xsn( 0 , 0 ) will 
diverge at a finite temperature for a finite and positive r independently how large the 
coupling constant g is. This means that a nematic transition will always take place 
before the long-range magnetic phase at r = 0 is reached if g is positive and T finite. 

In order to evaluate Eq. (J9]) we use the spectral representation for y, 

y(q ,iu n )=f de A ^ ^ , (11) 

J -oo - e 

with 

A(q, e) = - -- e -^— 2 -. (12) 

^7 [r + (q — Q q ) 2 ]^ + e 2 /q 2 

Inserting CD into (J9]) and carrying out the frequency sum yields after an analytic 
continuation for the imaginary part 

Im XsN(q> w ) = UJ / *A(k + q/2,e + o;)A(-k + q/2,e) [ra(e) - n(u + e)], (13) 
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with n(e) 
relation, 


l/(e e / T — 1). The real part of Xsn is computed via a Kramers-Kronig 


Rex^J(q,o;) = -p.v. 

7 r 


du 


Im xsN(q . v ) 

V — UJ 


(14) 


where p.v. denotes the principal value. The function xsN(q ! o’) is then obtained from 
(JHJ), (USD an d (ED - The spectral function of the spin nematic fluctuations is given by 


S{q.,w) = - [1 + n(w)]ImxsN(q,w), 

7r 

which we compute numerically. 


(15) 


3. Results 

3.1. Choice of parameters 

To compute the spin nematic spectrum numerically from (JHJ) and (T9|l . we first fix the 
spin fluctuation propagator [see (J6])] and the coupling strength g [see ©]• Neutron 
scattering data [30) for BaFei. 85 Coo.i 5 As 2 yield parameters for the spin fluctuation 
spectrum, namely, c ~ 1.3 and 7 ~ 230 meV. Hence we take in the present theory 

c = 1 and 7 = 1 , (16) 

measuring all quantities with the dimension of energy in units of 7 . For the mass term 
r we take in the normal state 

r = r cr + (T - T sn ) for T > T SN - (17) 

The bare susceptibility Xsn( 0,0) depends on T and r. The instability equation of the 
normal state with respect to the nematic state [see (I5T> ] . 

l = 0Xg$(O,O), (18) 

yields a relation between the critical mass r cr = t(Tsn) and the nematic transition 
temperature T SN . The value of T s n (or equivalently r cr ) can be considered as a free 
parameter and we will study the spin nematic spectrum for various choices of Tsn- At 
T = 0, QT7P is not valid and r should be taken as a non-thermal control parameter, for 
example, concentration of an isovalent substitution, carrier density, or pressure. 

According to the phase diagram of Ba(Fei_ x Co ;c ) 2 As 2 in [31], a structural phase 
transition and the SDW instability occur at 80 K and 70 K, respectively, at x — 0.04. 
If the structural phase transition is assumed to be associated with the spin nematic 
instability, the present theory reproduces the experimental transition temperatures at 
x = 0.04 for g = 0.29, extrapolating r linearly down to Tsdw- Since a different value 
of g would be obtained if one considers experimental data at a different x, we will also 
take g = 0.9 later to clarify how our results depend on the choice of g. 

We first study the case of g — 0.29, and calculate the critical value r cr as a function 
of T S n (figure [ 2 D and the spectral weight of spin nematic fluctuations in q and oj space 
(figures [3) HJ and ED- We will then consider results for a larger coupling strength g = 0.9 
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(figures EJ El and [7]) . On the basis of these results, we sketch typical phase diagrams of 
the spin nematic phase near the SDW phase (figure [8]) . 

3.2. Numerical results 

We put q — > 0 and oj —> 0 in (JHJ) and study the spin nematic instability. We search for 
its onset temperature Tsn by determining the temperature at which Xsn(0,0) diverges 
for a given r cr . Figure [2] shows r cr as a function of Tsn for two coupling strengths g. r cr is 
always positive and increases with increasing Tsn- This feature can easily be understood 
by using the approximate expression (fill for which holds at large temperatures. 
The instability condition Eq. (TT8|) becomes then equivalent to r cr ~ gfl sn- The linear 
dependence of r cr with Tsn is at least approximately reflected in the calculated curves 
in figure [2l The same holds for the increase of the slope with increasing g. 



Figure 2. (Color online) Critical value r CI as a function of the onset temperature of 
the spin nematic instability Tsn for a spin nematic coupling strength of g = 0.29 and 
0.9. 

Figure [3] shows the spectral weight of the spin nematic fluctuations at q = 0 and 
oj = 0 in the plane of Tsn and T—Tsn- The spectral weight is enhanced upon approaching 
Tsn and eventually diverges at Tsn- In particular, strong fluctuations appear at higher 
temperatures above Tsn if Tsn becomes larger, as seen in the yellow region in figure [3] 
For Tsn ~ 0, on the other hand, the enhancement of the spectral weight occurs only 
very close to T = T S n- 

Retaining still q = 0, we next present in figure |4] the dependence of the spectral 
weight on oj for several temperatures. At high temperature well above Tsn the spectrum 
is almost flat at low energies and its weight is very small. With decreasing temperature 
the low-energy spectral weight is enhanced to form a peak at zero energy in form of a 
central peak. This peak grows more and more upon approaching Tsn and finally diverges 
at Tsn- 

Figure [5] is a q-w map of the spectral weight near the spin nematic instability. The 
highest spectral weight is located around q = 0 and oj = 0 and the weight spreads 
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Figure 3. (Color online) Spectral weight of the spin nematic fluctuations at q = 0 
and lj = 0 for g = 0.29 in the plane of Tsn and T — Tsn; the spectral weight is plotted 
on a logarithmic scale. 



Figure 4. (Color online) u> dependence of the spectral weight at q = 0 for g = 0.29 
and different temperatures. The spin nematic instability occurs at Tsn = 0.05968 with 
the critical value r cr = 0.00968. 


with increasing energy while its strength is decreasing, like the tail of a comet. That 
is, the spin nematic fluctuations appear as a diffusive peak around q = 0 and u — 0, 
and no dispersive features can be seen. While we have chosen T S n = 0.05968 and 
T — Tsn = 0.01 in figure [5J. essentially the same result is obtained for other choices of 
parameters, although the tail of the comet is gradually blurred with increasing T. 

We next comment on results at Tsn = 0. Figure 2 shows that for the coupling 
strength g = 0.29 r cr almost vanishes at Tsn = 0. Hence the spin nematic and SDW 
instability occur almost simultaneously. Collective effects of spin nematic fluctuations 
occur only in the vicinity of T = Tsn as can be inferred from figure [3] As a result spin 
nematic fluctuations are enhanced only well below T & 10 -4 for Tsn — 0. If a larger 
value for the coupling strength g is used, spin nematic fluctuations for Tsn = 0 become 
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Figure 5. (Color online) q-w map of the spectral weight on a logarithmic scale for 
g = 0.29. The scanned direction of q is sketched inside the figure: (0.487T, 0.487r) 
—> (0,0) —> (0.487T, 0) —> (0.487T, 0.487t). The spectrum has a fourfold symmetry in 
momentum space above Tsn- 


visible at much higher temperatures. For g = 0.9, for instance, the overall temperature 
scale associated with the spin nematic instability becomes much larger than for g = 0.29 
and low-energy fluctuations are strongly enhanced already below T 0.02 [« 50 K, see 
(fT6l) and also figure [7]. We have checked also that there is no qualitative change in 
figures El El and [5] if a finite T SN is considered and g is increased from 0.29 to 0.9. 

At zero temperature we consider r in ([6]) as a non-thermal control parameter, see 
also the statement below (fl8|i . We plot in figure [6] the c o dependence of the spin nematic 
spectral weight at q = 0 for several values of r above the spin nematic instability at 
r cr = 0.00404 for g = 0.9. In contrast to the case of a finite T, described in figure El 
no central peak is formed and the weight at oo = 0 remains zero when approaching 
the spin nematic instability. Instead spin nematic fluctuations appear as a soft mode. 
With decreasing r a peak structure forms at a finite energy and moves towards lower 
frequencies. When r approaches r cr , the width of the peak becomes very narrow and at 
the same time the height of the peak increases strongly. At the spin nematic transition 
r = r cr the weight diverges at cu = 0. The presence of the soft mode suggests that there 
might be a dispersive feature of the spin nematic fluctuations in the plane of q and u at 
zero temperature. Hence we computed maps of the spectral weight, similar to figure El 
at T = 0. The dispersive feature is actually obtained, but only for r very close to r cr ; 
furthermore it becomes visible only in the vicinity of q = 0 and u = 0, on a much 
smaller scale than that in figure El In this sense, the dispersive feature is extremely 
weak. In fact, once r becomes a little larger than r cr , spin nematic fluctuations produce 
a diffusive signal around q = 0 and u> = 0 even at T = 0, similar as in figure El 
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Figure 6. (Color online) oj dependence of the spin nematic spectral weight at q = 0 
for several values of r. The weights for r = 0.0045 and 0.0042 are actually larger by a 
factor 2 and 7.14, respectively, as indicated in the figure. 


3.3. Asymptotic behavior of Xsn at low frequencies 

Our obtained results can easily be understood by analyzing the low-energy property 
of the spin-fluctuation bubble Xg^(q, u). From (TT3]1 follows Imyg^(q, — ui) = 
— Imxg^(q, uj), i.e., Inryg^q, uj) is an odd function of the frequency. Evaluating (1T3l) 
for uj, T -C r yields the asymptotic expansion 

rp2 

ImXsN(0> w ) ^ + «3 W 3 H-, (19) 

where <23 is a constant. There is a term linear in uj, which yields a central peak as found 
in figure U Its coefficient has a T 2 dependence, leading in figure [3] to strong fluctuations 
over a wide temperature region above T SN for a higher T SN . On the other hand, at T = 0, 
the liner term vanishes and the spectral weight is characterized by uo 3 at low uj, leading 
to a substantial suppression of the spin nematic fluctuations at low uj. This is the reason 
why spin nematic fluctuations are strongly suppressed at low temperatures (figure [3]) and 
no central peak is present. Instead a soft mode associated with spin nematic fluctuations 
occurs at T = 0 as shown in figure |HJ Since the spin-fluctuation propagator is in general 
parameterized by © close to the SDW instability, the low-energy dependence of (jT9|) is 
understood as a general feature of the spin-fluctuation bubble. This low-energy property 
yields characteristic features of the spin nematic fluctuations found in figures [3] - [ 6 ] as 
well as the spin nematic instability close to the SDW phase shown in figure [21 

The peculiar behavior of Xg^(0,u;) can be seen most clearly in figure [7] The two 
diagrams in this figure show the evolution of the spectral weight for g = 0.9 as a function 
of temperature. Tsn is equal to 0.00016 and thus very small, yielding spectra which are 
also representative for T SN = 0. Fixing T SN means that also r cr is fixed, see (Eli and 
figure El The lower diagram shows that for T > 0.02 the spectrum is completely flat 
and structureless. Decreasing T down to 0.004 the spectral weight shifts towards lower 
frequencies and a central peak is formed. In this temperature range the linear term 
in uj of lmxg^(0, of) still dominates in a low frequency expansion of Imy^ [see (1X9]) ] 
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Figure 7. (Color online) Evolution of the spectral weight for g = 0.9 and Xsn = 
0.00016 as a function of temperature. The critical value r CI is 0.00406 (see figured]). 


and causes the central peak. Lowering further T to 0.001 the low frequency part of 
the spectral weight looses intensity and a propagating peak appears with an increasing 
frequency and a decreasing half-width. In this temperature interval the linear term in 
c o of lmxg^(0, a>), proportional to T 2 , becomes small and the term a^u: 3 starts to play 
a role. As a result a transition from a diffusive to a propagating mode behavior is 
obtained. Considering the upper diagram and decreasing further the temperature the 
energy of the propagating peak and its half-width decrease whereas its height increases 
strongly. 

Sfl. Schematic phase diagram 

On the basis of our obtained results (figures [2] and [3]), we can infer a typical phase 
diagram for the spin nematic and SDW phases. Let the onset temperature of the SDW 
instability evolve as a function of a control parameter 6, as shown in figure [8](a); 6 
may correspond to the concentration of substituted ions, carrier density, pressure, or 
other quantities depending on material properties. Xsdw decreases with increasing 5 and 
vanishes at a critical value of (>sdw- Let us introduce the slope a = r cr /(T SN — T SDW ) 
(> 0) and assume that a depends only weakly on 5. Approximating the results in 
figure [2] by r cr oc T sN , we obtain T sN — T SDW oc T SDW , i.e., the temperature region 
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occupied by the spin nematic phase increases with increasing T s DW . In the opposite 
limit where T s D w is close to zero, i.e., 5 —> <5sdw; the spin nematic phase vanishes very 
close to $sdw- As seen from figure [3l the temperature region where strong spin nematic 
fluctuations are present expands for a higher Tg N and shrinks substantially near <5sdw as 
shown by the dashed line in figure El^a). On the other hand, if the strength of the spin 
nematic interaction becomes very strong, the spin nematic phase as well as the region 
where strong spin nematic fluctuations are present expands to a larger region as shown 
in figure E](b). In particular, a region of the spin nematic phase can be well separated 
from the SDW instability at zero temperature (see also figure [2] for g = 0.9) and thus 
two well-separated quantum phase transitions exist as a function of 5. 



Figure 8. Schematic phase diagram of the spin nematic (SN) instability near the 
SDW phase on the plane of a control parameter S and the temperature for a realistic 
value of g for iron pnictides (left figure) and for a large g (right figure). The phase 
transition is assumed to be continuous. Two quantum phase transition occur at <5sn 
and ^sdWj and these two almost coincide for a small g. 


4. Discussions 

4-1. Origin of the nematic phase in iron-based superconductors 

We discuss the origin of the nematic phase observed in iron-based superconductors (8j. 
In figure 0 we assumed that Tsdw decreases monotonically with increasing S, which is 
the usual case in iron-based superconductors. The spin nematic scenario then predicts 
that the temperature difference of T SN and T SDW , namely r cr /a — Tsn ~ T sdw > 0, 
should decrease monotonically with increasing <5, at least, if a depends only weakly on 
6. However, this tendency is not observed in iron-based superconductors in spite of 
the fact that many different compounds [7] have been investigated. We discuss several 
possibilities to resolve this qualitative discrepancy between the expectation and the 
experiment. 
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First, we have assumed a constant coupling strength g and a constant value of c; 
the latter controls the overall strength of spin fluctuations [see (jBJ)]. If g (and/or c) 
is substantially suppressed at a low <5, Tsn would shift closer to Tsdw similar to the 
experimental observation. However, a theoretical study [29] suggests that the value of 
g becomes larger with decreasing <5, so that r cr is expected to be larger than in figure [8] 
at low 6. The discrepancy between the present theory and experimental observations 
would even increase. In addition, the value of c is expected to become larger for a lower 
6 because the nesting condition of the Fermi surfaces for the momenta (tt, 0) and (0,7r) 
becomes better at a lower 5. Consequently, a more realistic treatment of g and c would 
lead to a larger discrepancy between theory and experiment. 

Below Tsn the spin nematic order parameter <j>(T)(> 0) becomes nonzero. As a 
result r(T) has the form ± 0(T) where r^°^(T) is equal to r cr + (T — Tsn) in 

agreement with the expression (fTTl) in the normal state. The temperature dependence 
of r will approximately be given by the solid curve in figure [9] The sudden drop 
of r just below Tsn is due to the development of the spin nematic order parameter 
characterized by <f>(T) (TgN — T) 13 ] (3 = 1/2 in mean-field theory and (3 = 1/8 in a 
two-dimensional system. The red line in figure [9] is a linear interpolation through the 
points (T, r) = (Tg N , r cr ) and (T SDW , 0) with the slope a. The absolute value of a simply 
changes the scale of temperature and thus is not relevant to the present discussion as 
long as a is independent of 5. 



Figure 9. (Color online) Temperature dependence of r = (£ is the magnetic 

correlation length) in the actual system (solid curve) and in the linear approximation 
(red line). The values of Tsn, Tsdw and r cr are correctly captured in the linear 
approximation, although the whole temperature dependence of r is not. The slope of 
the red line corresponds to a and is assumed to depend weakly on 5 in figure [8] 

To reconcile the observed nematic phase in terms of the spin nematic instability, we 
have to invoke a 5 dependence of a. In the so-called ”1111” compounds LaFeAsOi_ x F x 
[32] and CeFeAsOi-zF^ [33] the onset temperature of the nematic phase largely extends 
to the superconducting region even though the SDW state has already vanished. A 
similar feature is observed also for Ba(Fei_ a ,Co a; ) 2 As 2 [31], where the nematic phase is 
confined closer to the SDW phase as compared to the ”1111” compounds. To understand 
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these data, a must be assumed to become very large near x ~ 0 so that T SDW occurs 
much closer to T SN . It also should decrease substantially with increasing x so that T s D w 
occurs further away from Tsn- This S dependence should be so dramatic that it fully 
changes the qualitative features of the phase diagram (figure [8]) and also compensate the 
possible further discrepancy caused by a realistic treatment of g and c. While the actual 
SDW transition changes to first order close to x = 0 whereas the nematic transition is 
continuous there [34], such a region is rather small and thus does not modify our major 
discussion. In typical hole-doped compounds Bai_ ;E K x Fe 2 As 2 [351 EE] and isovalent 
doping materials Ba(Fei_ 3 ,Ru a .) 2 As 2 [371 38] the nematic phase occurs simultaneously 
with the SDW instability within experimental resolutions. To understand this, a should 
be assumed to be very large near x = 0 and to remain rather large with increasing x. 

At present it is not clear whether the assumption of a 5 dependence of a is justified 
from a microscopic point of view. To evaluate a the absolute value of 0(T) does matter 
as a function of T and moreover an approximation scheme to calculate r also does 
matter since r becomes zero only at T = 0 in a purely two-dimensional system. Because 
of these subtleties in determining the precise form of the spin nematic phase, it seems 
natural to expect a substantial material dependence of the shape of the spin nematic 
region in the 5-T plane. Thus one would naturally expect a phase diagram similar to 
figure [8] at least for a certain class of materials if the spin nematicity is responsible for the 
nematic phase in iron-based superconductors in general. At present, however, no iron- 
based superconductors are known to show characteristic features of the spin nematicity 
shown in figure E] [ZJ. Fernandes et al. study a possible phase diagram of the spin 
nematic phase near the SDW phase [29]. Their calculations are done in two limits: the 
zero-temperature limit and the classical limit in the sense that only the zero Matsubara 
frequency is considered in gap equations. Their results at T = 0 predict that the spin 
nematic and SDW instabilities occur simultaneously. This is consistent with our results, 
although they predict a first order transition, a possibility which is not considered in 
our analysis. At finite temperatures they obtain in the classical limit figures 8 and 14 in 
[29] which compare successfully with experimental data and also present a microscopic 
treatment of the slope a. In our approach all Matsubara frequencies are kept but a is 
considered as a phenomenological input. 

Our schematic phase diagram (figure [8]) does not apply to the orbital nematic 
scenario, because the orbital nematic instability is controlled by a fermionic bubble 
diagram [391 Mi- Hence it is tempting to state that in general orbital nematicity is 
likely responsible for the nematic phase observed in iron-based superconductors [El IT]- 
However, it seems too early to reach such a conclusion. In view of the fact that the 
nematic phase occurs in general close to the SDW phase, it is important to clarify from 
a microscopic point of view why this should also hold for the orbital nematic state. 
One possible reason is given in [5] where Kontani et al. point out the important role of 
Aslamazov-Larkin type diagrams. It is interesting to explore further whether the orbital 
nematic scenario indeed explains the fact that the nematic instability occurs close to the 
SDW phase at 6 = 0 and the nematic region extends to a larger region with decreasing 
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Tgow and higher 8 as observed in experiments HI Eli- 

The iron-based superconductor FeSe shows a structural phase transition from an 
orthorhombic to a tetragonal phase with decreasing temperature, but no SDW phase has 
been detected [30]. This experimental fact naturally suggests that orbital nematicity is 
responsible for the structural phase transition. If we wish to understand the structural 
phase transition in FeSe in terms of the spin nematic order, we would invoke, for example, 
a large coupling strength g, as shown in the right panel in figure [8] FeSe would then be 
located in the region 5 sdw < 8 < 5sn- 

4-2. Nematic fluctuations 

Spin nematic fluctuations can be measured directly by electronic Raman scattering. 
Although the computation of the Raman intensity in a microscopic model for electrons 
is beyond the scope of the present study, our obtained spectra can be interpreted as 
B\g Raman spectra within the following approximation. In a microscopic calculation 
the vertex 7 k in figure DJa) is replaced by a triangle diagram constructed from three 
electronic Green’s functions. Such a triangle diagram depends both on momentum 
and frequency. Expanding its momentum dependence in terms of a complete set of 
functions with B\ g symmetry it is evident that only the function y(k) contributes 
substantially due to the restriction of the momenta to the neighborhood of (n, 0) 
and (0,7r). Concerning the frequency it is plausible that the triangle diagram has no 
resonances at low frequencies in the range of collective nematic fluctuations. Thus the 
triangle diagram may be approximated by a constant and we expect that our results 
in figures [MZI capture the major features of electronic Raman scattering due to spin 
nematic fluctuations in the Bi g channel. Although available Raman scattering data 
mmm are interpreted in terms of orbital nematic fluctuations |39], they do not 
seem to exclude the spin nematic scenario. Not only further theoretical studies but also 
more detailed experimental data are important to determine the origin of the nematic 
phase observed in iron-based superconductors. A crucial test to identify the origin of the 
nematic phase is to measure nematic fluctuations at zero temperature by suppressing 
the superconductivity, for example, by applying a large magnetic field: spin nematic 
fluctuations appear as a soft mode upon approaching the nematic phase whereas orbital 
fluctuations form a central peak. 

We finally discuss a possible role of nematic fluctuations for superconductivity. If 
the nematic phase in iron-based superconductors is associated with spin nematicity, 
one might wonder about superconductivity mediated by spin nematic fluctuations. As 
shown in figures [3] and [8l spin nematic fluctuations are substantially suppressed at low 
temperatures. It is thus unlikely that such fluctuations are responsible for the observed 
superconductivity. Instead usual spin fluctuations provide a more natural scenario to 
understand superconductivity mm in these sytems even if the spin nematic phase 
occurs in actual materials. On the other hand, if the nematic phase originates from 
orbital degrees of freedom, it has been shown that orbital nematic fluctuations can 
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lead to strong coupling superconductivity [2HJ in pnictides and thus provide an exotic 
mechanism for superconductivity in these systems. 

5. Conclusions 

Using a general form for the spin-fluctuation spectrum near the SDW phase we 
have studied spin nematic spectra in energy and momentum space. If the spin 
nematic interaction is attractive and the nematic transition continuous, its transition 
temperature Tsn is always higher than that to the antiferromagnetic state at Tsdw- 
The spin nematic spectra are characterized by several generic features: (i) the critical 
magnetic correlation length £ cr decreases with increasing Tsn (figure [2]), (ii) strong low- 
energy spin nematic fluctuations extend to a wider temperature region for a larger 
Tsn (figure [3]), (iii) the spin-nematic spectrum at q = 0 exhibits a central peak as a 
function of to close to the spin nematic instability at a finite temperature (figures 0] and 
ED whereas it shows a soft mode upon approaching the spin nematic instability close to 
zero temperature (figures [6] and ED , and (iv) there is no clear dispersive mode associated 
with spin nematic fluctuations and instead a diffusive peak is obtained around q = 0 
and uj — 0 (figure ED . These general features originate from the low-energy property of 
the simple bubble diagram of spin fluctuations [see (TT9l) ] . The resulting phase diagram 
is shown in figure [8] if a , the average slope of r cr = £“ 2 as a function of temperature 
between Tsdw and Tsn, can be considered as a constant. The existing discrepancies 
with the experimental phase diagram may indicate that a constant a is not adequate 
or that another scenario such as orbital nematicity may be more appropriate. Since a 
rather general form for the spin-fluctuation spectra [see (J6J)] is employed in our study, 
the present theory can be applied or extended straightforwardly to other systems where 
magnetic fluctuations are characterized by four wavevectors with fourfold symmetry. 
Although cuprate superconductors exhibit nematicity in the magnetic excitation spectra 
[44] and thus might be possible systems for spin nematic order, the line of onset 
temperatures of nematicity is nearly parallel to that of the incommensurate magnetic 
order as a function of doping in YE^CusOe+y [45], which is at variance with figure IH1 
Instead, the nematicity in cuprates was discussed in terms of a feedback effect from 
charge nematicity [06J 0?] . 
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